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1 Introduction 

Numerous studies of the somatosensory evoked field 
(SEF) have confirmed the homuncular organization 
of the first somatosensory cortex (SI) [1, 2]. A 
dipole source analysis of SEFs recorded near the 
surface of the head allows to localize the underlying 
neural generators in SI and in the secondary 
somatosensory cortex (SII) [3, 4]. Theoretical as 
well as experimental studies suggest that sources can 
be localized relative to each other with an accuracy 
of the order of millimeters [5, 6]. In the case of the 
auditory evoked field (AEF) single time slices can 
generally be described quite well by an equivalent 
current dipole (ECD). However, the locations 
derived in this way generally exhibit a considerable 
movement along a trajectory [7]. Purpose of the 
present analysis was to find out whether a similar 
spatio-temporal structure of the underlying sources 
can also be found in SI. 

2 Methods 

2.1 Subjects 

Neuromagnetic measurements were done with a 37 
channel first-order gradiometer system (Magnes, 
Biomagnetic Technologies Inc., San Diego, USA) in 
a magnetically shielded room. The subjects were 
comfortably lying on their right side. The magnetic 
field was recorded over the contralateral hemisphere 
with the gradiometer array centered between the 
positions of C3 and T3 according to the standard 
10/20 system (see Fig. 1). 



Figure 1: The positioning of the gradiometer system. 


2.2 Stimulation and experimental procedure 

SEFs were measured by applying tactile stimuli to 
the digits 2 to 5 (D2, D3, D4, D5) of the right hand 
using balloon diaphragms driven by bursts of 
compressed air (see Fig. 2). The stimulus duration 
was about 150 ms. The stimuli were presented in 
randomized order with an interstimulus interval 
(ISI) of 0.5 s and 1 s, respectively [8]. Between two 
stimulations of the same digit at least one 
stimulation of a different digit was interspersed. 
Thus, the minimal interval between two stimulations 
of the same digit was twice the ISI, whereas the 
respective mean interval was four times the ISI. The 
total number of stimuli was basically the same for 
every digit. Four independent measurements were 
performed for both ISIs [9]. 

The subjects were asked to maintain their position 
throughout the runs and to pay no attention to the 
stimuli. To facilitate the observance of the latter 
point, they were watching a video without sound. A 
total number of eight right-handed subjects (age: 24- 
32 years) with no history of neurological or 
psychiatric disorders was investigated. 



Figure 2: This figure shows the balloon diaphragms 
placed on the finger tips. During the actual 
measurements the hand was kept more relaxed, 
though. 

2.3 Data analysis 

The different measurement conditions were 
compared by inspecting the root-mean-square 
(RMS) values of the measured fields as well as the 







moments derived by modeling the data with a 
moving ECD in a sphere. 

In the time interval from 30-60 ms after stimulus 
onset only those sources have been taken into 
account for further analyses which were able to 
explain the SEFs with a goodness of fit (GoF) value 
greater than 90%. 

3 Results 

In all subjects, the estimated moving dipole 
locations described reproducible trajectories, 
generally showing little overlap when considered in 
all three dimensions. Their length was typically of 
the order of 5 mm which was to some extent 
dependent on the GoF criterion used (see Fig. 3). 

In agreement with an earlier study, the IS I had no 
significant influence on the source locations [9]. 



Figure 3: Source trajectories for digits 2 to 5 in the 
latency range from 30-60 ms after stimulus onset 
are shown in the x-z-plane. 


The geometrical arrangement of the trajectories 
showed a great interindividual variability. However, 
in most subjects they appeared to run more 
perpendicular rather than parallel to both the 
orientation of the dipole moment and the 
somatotopy axis, which was defined as the line 
connecting the source locations with the maximum 
RMS value for each digit. 

All sources of the trajectories had in common that 
they were located in the SI cortex which could be 
verified by overlaying the estimated dipole locations 
with magnetic resonance images (see Fig. 4). An 
exemplary 3-dimensional view of the trajectories is 
shown in Fig. 5. 



Figure 4: The three figures above show overlays of 
several MR images per view to display all source 
locations of the trajectories which penetrate several 
MRI slices. The sources are displayed as small 
overlapping circles which are all located in SI. 






Figure 5: Exemplary 3D-plot of source trajectories 
of a single run. The four digits are color coded as 
follows: 

D2 = red, D3 - green, D4 - cyan, D5 - blue. 

The coordinate system indicates the orientation of 
the figure: x points to the right, y up and z to the 
front. The yellow line represents the somatotopy 
axis. 

The direction of propagation of the moving dipoles, 
which are describing a single trajectory, was defined 
as the vector from a dipole location about three time 
steps before the latency of the dipole with the 
maximum RMS value (D max ) to another dipole 
location about three time steps later. Calculating the 
scalar product of this direction with the orientation 
of D max revealed that the dipole orientation is almost 
perpendicular to the direction of the dipole 
movement (see Fig. 6). 
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Figure 6: The histogram shows the distribution of 
angles between the orientation of D max and the 
somatotopy axis (red) as well as the angle between 
the orientation of D max and the direction of source 
movement (blue). 


Furthermore, the scalar product of the direction of 
the dipole movement with the somatotopy axis 
showed that the angle between them varies mainly 
around 70°-85° degrees. 

The distances D2-D3, D3-D4, D4-D5, and D2-D5 
between the activated cortical areas were calculated 
for each digit from the locations of D max . The 
median values for these distances averaged over all 
subjects are listed in Table 1. 


Table 1: Median distances between cortical areas 
activated by stimulation of digits 2 to 5. 


digits 

median distance 

D2-D3 

0.39 cm 

D3-D4 

0.16 cm 

D4-D5 

0.38 cm 

D2-D5 

0.48 cm 


4 Discussion and conclusions 

Compared to the results obtained from ECD 
calculations representing the source of the auditory 
Nlm component [7], the ECD representing the SI is 
relatively stable. Nevertheless, trajectories were 
found also in the present study. The fact that the 
estimated dipoles (assumed to have an orientation 
perpendicular to the cortical surface) are roughly 
perpendicular to these trajectories, suggests that the 
trajectories describe spreading activation along the 
cortical surface. The large interindividual variability 
of the dipole trajectories reflects the anatomical 
variability of the postcentral sulcus. 

Since the somatotopy axis was defined by a line 
connecting the locations of D max for each digit, it is 
not a straight line, but a piecewise linear one which 
approximately follows the cortical surface. 
Therefore, the angles between this axis and the 
orientation of the dipoles D max , presumably 
perpendicular to the cortical surface, have a great 
variability and are in most cases slightly less than 
90° degrees. An additional point is that only the 
tangential component of the magnetic field has been 
taken into account for ECD calculations since the 
MEG is insensitive to radial components. 

The sum of the Euclidean distances D2-D3, D3-D4, 
and D4-D5 turned out to be about twice the length 
of D2-D5. Thus, the latter measure, which does not 
account for the folding of the cortex, characterizes 
the extent of the SI representation of the digits only 
insufficiently. 
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